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Abstract 

The flow of amorphous solids results from a combination of elastic deformation and local struc- 
tural rearrangements, which induce non-local elastic deformations. These elements are incorporated 
into a mechanically-consistent mesoscopic model of interacting elastoplastic blocks. We investigate 
the specific case of channel fiow with numerical simulations, paying particular attention to situ- 
ations of strong confinement. We find that the simple picture of plastic events embedded in an 
elastic matrix successfully accounts for manifestations of spatial cooperativity. Shear rate fiuctua- 
tions are observed in seemingly quiescent regions, and the velocity profiles in confined flows at high 
applied pressure deviate from those expected in the absence of nonlocal effects, in agreement with 
experimental data. However, we suggest a different physical origin for the large deviations observed 
when walls have rough surfaces, associated with "bumps" of the particles against the asperities of 
the walls. 
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Shear waves are damped in liquids, and propagate in elastic solids. This well-established 
distinction has major implications in the field of seismology but also has bearing on the 
peculiar nonlocal rheology of soft jammed/glassy materials, which share solid (elastic) and 
liquid (fiow) properties. 

In fact, the fiow of these materials bears notable similarities with earthquakes: it features 
a solid-like behavior at rest and local yielding above a given applied stress. Yielding is 
characterized by the emergence of local "shear transformations" [2] involving a few particles 
associated with a local fiuidization of the material. These structural rearrangements, 
hereafter named plastic events, induce long-range deformations. The microscopic details 
vary to some extent with the particular nature of the material. In the case of foams, they 
are identified as Tl events [4] in which the local change of neighbors is mediated by an 
unstable stage with four bubbles sharing one vertex. The robustness of the scenario for an 
extremely wide range of materials is striking. Ample evidence of the local plastic events 
and their long-range eflFects is indeed provided both by experiments using diverse materials 
0, y, gI and by simulations 0, Isl . 

Turning to the specific case of channel fiow, the existence of long-range interactions render 
any purely local approach questionable in a system featuring an inherently inhomogeneous 
shear and specific boundary conditions due to the walls. Convincing experimental evidence 
disproves the very existence of a local constitutive relation a = a (j) [9], where 7 and a are 
the local shear rate and shear stress. The presence of nonlocal eflFects highlights the eflFect 
of spatial cooperativity. These observations were rationalized in j^, [lo| by introducing a 
diflFusion equation for the local fiuidity f = ^/a: 

eM-{f-huik{<y)) = ^, (1) 

Here is a cooperativity length, and fijuik is the local bulk fiuidity, measured in a homo- 
geneous, simple shear situation. EqU] has been tested with considerable success 

HQ. 

However, the approach ignores the fiuctuating nature of plastic deformation, and relies on 
experimental measurements of the boundary conditions. 



Mesoscopic models |l3l-ll5| oflFer a computationally-eflBcient way to recover part of the 
complexity of the dynamics, while leaving behind microscopic details. Most importantly, 
they allow us to test our understanding of the physical processes involved in the fiow, and 



identify relevant parameters for its description [jj 



. However, direct comparison with ex- 
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perimental data is scarce. 

In this Letter, we build upon previous work [l^,[l5| to develop a simple, but mechanically 
consistent, 2D tensorial mesoscopic model that incorporates the phenomenology reported 
above, as well as (coarsened) convection. For the first time such a model is used to simulate 
a channel fiow, and its predictions are compared to experimental data. 

The channel is modeled as a rectangular elastic matrix, which is spatially discretized into 
a regular lattice of square-shaped blocks. Under the assumption of isotropy and incompress- 
ibility, Hooke's law states that the deviatoric stress tensor cr is proportional to the local 
strain tensor e, 

^xx M i^^xx ^yy) ^l-^^xx ) ^xy ^l-^^xy (2) 

/i is the shear modulus, and x and y are the fiow and gradient directions, respectively. The 
system covers the domain (x,^) G [0, L^^] x [Q^Ly]^ with and Ly the channel length and 
width. The initial response of the material confined in a channel to an applied pressure 
gradient Vpe^ is therefore given by: a^y (x, y^t = 0) = Vp {y — and Gxx (^5 V^t = 0) = 
0. 

As soon as a yielding criterion is met in a block, the block undergoes a plastic event, 
during which the stored elastic energy is dissipated. Within the framework of eigenstrain 
theory [16], the elastic strain e is converted into a plastic eigenstrain e^^ Physically, this 
corresponds to the change of reference elastic configuration following particle rearrangement. 
The conversion takes a finite time, due to dissipative effects. 

Since the block is embedded in an elastic matrix, the localized eigenstrain e^^ induces an 
elastic field cr^-^^ (^) = f G (^5 ^0 ' (r^) dr\ where the elastic propagator Q depends on 

the boundary conditions. The assumption of an underlying elastic behavior is bolstered by 
the similarity of the induced field with experimental and numerical observations jsl, [sl , all of 
which exhibit a quadrupolar angular dependence and a ||r — r'||~^-like decay in d dimensions. 
Combining the different elements, one gets the equation of evolution of a block: 

dta (r) = S^^* {r) + J G (r, • (/) rfV' (3) 

Here, e^^ = if the block is plastic, e^^ = otherwise, and r is the characteristic timescale 
for the release of the stored elastic stress. S^** is the time-derivative of the stress response 
of a purely elastic material to the same applied conditions; therefore S^^* = if the fiow is 



pressure- driven. Numerically, at each time step, Eql3] is solved in Fourier space, on a finer 
grid where, for accuracy, each block is divided into four subcells. 

There is obviously some arbitrariness regarding the choice of the yielding criterion and 
the duration of plastic events. The corresponding switching rates are denoted by / (a) and 



The probability to yield / (a) is set to / (a) = B (a — a^y) exp y-^j where B is the 
Heaviside function and a = ||cr||, so as to mimic an activated process of non-cooperative 
origin, associated with a material-dependent intensive parameter xi. Imposing a finite crit- 
ical stress cr^^, below which no plastic event can occur, is crucial for the existence of a 
macroscopic yield stress a (7 ^ 0) 7^ 0. Also note that the von Mises yield criterion is 
recovered in the limit xi ^ 0. The above expression accounts for the long-time relaxation 
of a pre-sheared material, observed even in the case of (athermal) granular materials |l7|. 
We assume that elasticity is restored when the plastic deformation rate becomes low, and 
choose: e (a) = exp where we have introduced another material parameter, Xg. 

Walls are modeled by imposing no-slip boundary conditions at ^ = and y = Ly. This 
results in a lengthy expression for the propagator, which we defer to a longer publication, 
along with details of the derivation. It is however noteworthy that, for a given eigenstrain 
the elastic stress that is released locally is up to 35% larger than in the bulk if the plastic 
event occurs near a wall. 

Finally, a coarsened picture of convection is implemented, whereby lines of blocks in the 
fiow direction (instead of individual blocks) are incrementally shifted. 

We choose units of time and stress such that r = 1, /x = 1 and we set ay = 1. The 
remaining model parameters, a^^, X/, and Xg, are now fitted by comparing the predictions 
of the model in simple shear fiow, i.e., a (t = 0) = and S^^^ = /n^app in Eql3l where ^app 
is the applied shear rate, with the macroscopic rheology measurements collected by Goyon 
and co-workers jsl, Q for an oil-in- water emulsion of average radius 6.5 /xm , which follows a 
Herschel-Bulkley law a = a^ + A^y^^ n ^ 0.5. The inset of Fig|3]shows that the experimental 
data can satisfactorily be reproduced with our model by using the following parameters: 
<^/iy=0.17, xi = 0.249, and x^ = 1.66. Note that the model units of stress and time have 
been appropriately rescaled. 



e(a) 



1(a) 

elastic regime ^ plastic event. 
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Cooper at ivity is a general feature of the flow of amorphous sohds, regardless of the flow 
geometry and the driving. Channel flow, however, is speciflc in that 

(i) the non-locality of the stress redistribution couples streamlines subject to different 
shear stresses. 

(ii) the presence of a wall, whether it be rough or smooth, may create a speciflc surface 
rheology, different from that in the bulk. 

In practice, these effects are of primary importance for conflned flows in microchannels, 
but they are not intrinsically caused by conflnement. 

Let us flrst consider Point (i). To estimate its importance, we introduce a dimensionless 
number, the Babel number Ba = S f / f^uiki^) ^ where Sf = f — f^uik is the deviation from 
the expected fluidity proflle owing to cooperative effects between regions subject to different 
driving forces. Within the theoretical framework of fluidity diffusion, Eq. [D, and for a fluid 
obeying a Herschel Bulkley relation in the bulk, it can be shown that Ba = ^(^| 



iivair^ 



• (J-(J0 



Large deviations from bulk behavior are expected when the stress is close to the yield 
value. A striking consequence of this flrst point was very recently unveiled by Jop et al. 



Uj in a system similar to that of Goyon et al. [l2j. They analyzed experimental data 
for an oil-in-water emulsion in an almost two-dimensional microchannel flow, and showed 
that the seemingly quiescent, rigid plug is in fact subject to flnite shear rate fluctuations 
^7 (x, y) = (7 (x, y)^) — (7 (x, y))^, where the brackets denote time averages. This clearly 
points to a nonlocal effect of the sheared regions. 

We now compare the model results[18] directly to experimental data: FigH] shows the 
velocity proflles [l^, and Figl2] shows the corresponding shear rate fluctuation proflles, for 
a numerical channel of transverse size A^^ = 16 blocks. 

Although the parameters of the model have been constrained only weakly by adjusting 
the macroscopic flow curve, semi-quantitative agreement is observed in regions far from the 
walls - apart from the large discrepancy at the highest applied pressure. The discrepancies 
in the vicinity of the walls will be considered below. It is interesting to note that the fltted 
channel size provides an estimate for the linear size N of an elastoplastic block measured 
in particle diameters, A^ 2. 

Not only shear rate fluctuations, but also the average shear rate, elude a local, or bulk, 
description, when Ba gets large. This paradigm-shifting fact is the major result of experi- 
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Figure 1. Offset velocity as a function of the channel crosswise coordinate, for stresses at the wall 
=141 Pa, 188 Pa, 235 Pa, 282 Pa, corresponding to = 0.36, 0.48^ 0.60, 0.72 in model units, 
from top to bottom, (x) Experimental data collected by Jop et al [Ul] (channel width: 225 /xm), 
(•) numerical results for a channel of Ny = 16 blocks crosswise. The curves are offset with respect 
to each other for clarity. 




Figure 2. Shear rate fluctuations (y) (averaged along the x-direction), for =141 Pa, 188 Pa, 
235 Pa, 282 Pa (identical to FigH]), from bottom to top. (x) Experimental data collected by Jop et 



al 



111, {solid lines) numerical results for Ny = 16. 



mental observations by Goyon and co-workers in microchannel flows of emulsions with either 
smooth or rough walls [9]. 

In the case of smooth walls, deviations from the macroscopic bulk constitutive relations 
are comparatively weak. In contrast to strain rate fluctuations, they become perceptible 
only at high applied pressures in confined geometries, when the stress at the wall is larger 
than some channel-size dependent critical value a^. Under these conditions, the overall 
fluidity of the system is enhanced. Consistently with our expression of the Babel number, 
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Figure 3. Velocity profiles in the micro channel. The curves are shifted for clarity, (x) Experimental 
measurements by Goyon et al. jl^l for (from top to bottom) (jy^ — 45, 60, 75, 91 Pa, i.e., cr^ = 
0.75, 1.0, 1.25, 1.52 in model units, in a 112/xm-wide smooth channel; (•) numerical results for 
Ny — 7. (The solid line is a guide to the eye); dashed black line: predicted velocity profile on the 

)asis of the macroscopic flow curve S {^app)^ shown in the inset graph ((x) experimental data from 

I2I 1: (•) numerical results in a large system, Ny — 64). 

the most prominent discrepancies with the bulk predictions are observed near the edges of 
the plug, which become more rounded. Accordingly, cooperativity can soften the material 
so that it flows even below the macroscopic yield stress. This aspect is corroborated by 
numerica. 
equation 

critical wall stresses to coincide between the model and the experiments, the transverse 
size of the numerical channel must be set between 6 and 10 blocks, which again corresponds 
to N0 ^ 2. Besides, for such transverse microchannel sizes, the maximal velocity displays 
large oscillations in time, as has often been reported experimentally [21]. 

Let us now turn to a comparison with a description in terms of fluidity diffusion, Eq. [H To 
solve Eq. [H the shear-rate dependence of the length ^ must be specified, and two boundary 
conditions (BC) are required. Regarding the BC, we impose f {y = 0) = f {y = Ly)^ and set 
the fluidity at a point close to the wall to the value measured in simulations. For the fluidity 



simulations [20I, and is expected within the framework of the fluidity diffusion 
lo| . Fig |3] shows that our model also satisfactorily captures these features. For the 



^, two cases are studied in FigHl either no dependence, i.e., C = Co, following 
1^ or a power-law dependence, ^ (7) = <^o7~°'^^5 where 7 is the product of the 



dependence o: 
Goyon et al. 

local shear stress and fluidity, as derived in Ref. [10] in the limit 7^0, and in reasonable 
agreement with the data of Ref. [11]. In both cases, ^0 is adjusted by a least square 
minimization. Overall, the power-law dependence provides a closer fit of our numerical 



0.12 
0.10 
0.08 
0.06 
0.04 
0.02 
0.00 

0.0 0.2 0.4 0.6 0.8 1.0 

Ly 

Figure 4. Fluidity profiles for Ny — 12. Filled circles: numerical results, dashed green line: solution 
of EqlDwitli i (7) = 0.03702, solid blue line: solution of EqlDwitli ^ (7) = 0.011467"°-2^ The thin 
dash-dotted lines represent the bulk fluidity fhulk- 

results, especially at higher applied pressures. However, neither assumption concerning 
^ (7) was able to provide a perfect fit, a defect that we ascribe to the approximation of 
long-range interactions by a diflFusive term, and to the neglect of fiuidity fiuctuations. 

Now, when rough walls are substituted for smooth walls, the situation diflFers widely. 
Much larger deviations are observed, even at lower Babel numbers than above. The discrep- 
ancies caused by the change of the surface roughness of the walls point to the prevalence the 
specific surface rheology, point (ii) above. Further evidence comes from the fact that devi- 
ations are particularly large close to the walls. The no-slip BC at the wall imposed in our 
model are insuflScient to account for these large deviations. One is therefore led to conclude 
that either (a) the roughness of the wall alters the structure of the material in its vicinity 
or (b) slip along the wall generates a stress field in the system, as particles constantly bump 
into wall asperities. 

Since the observed deviations are not restricted to the direct vicinity of the (rough) 
walls, we examine the second possibility here, namely slip along the wall as a source of 
mechanical noise. Goyon et al. [9] measured wall slip both for smooth and rough surfaces, 
the latter having asperities of a characteristic lengthscale of 1 /xm, whereas the emulsion is 
made of ^ 6 /xm-diameter droplets. In the case of smooth surfaces, asperities are too small 
to deform the particles significantly, and the presence of wall slip is therefore unlikely to 
induce considerable change as compared to the no-slip case, which validates our approach. 
On the contrary, wall slip along a rough surfaces brings about "bumps" (or deformations) 
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of the particles into surface asperities. Similarly to bulk plastic events, these "bumps" are 
expected to induce an elastic field in the material. This scenario has the potential to explain 
why deviations are not always observed in nominally similar conditions: Seth and co-workers 



2211, instance, used a different surface preparation protocol to obtain rough surfaces and 
did not detect any nonlinearity of the fiow in the vicinity of the rough walls, but, quite 
interestingly, they also reported that wall slip remained negligible in that case. 

Our model cannot render the mechanical effects of wall slip along a rough surface without 
further input. Nevertheless, fictitious plastic events can artificially be added along the walls 
to mimic the impact of "bumps". Concretely, we now model both walls as lines of blocks and 
select a fraction (1/3) of these blocks at random as mechanical noise sources, that is to say, 
they shall release a constant plastic strain e^^y^^^ {'^ 5) by unit time. It should be remarked 
that the procedure does on no account violate mechanical equilibrium. The resulting local 
flow curves, shown as Supplemental Material, are qualitatively similar to those obtained 
experimentally in Ref. (qI. Let us also remark that shear rate profiles in the presence of 
fictitious plastic events do not fiatten in the vicinity of the walls {data not shown) ^ as in 
Figl2l and are therefore more compatible with the experimental data. 

In conclusion, we have refined and extended a mesoscopic model based on a default 
elastic behavior and finite-time plastic events to simulate a confined channel fiow. We 
have realized the first direct comparison of such a model with experimental data, paying 
particular attention to manifestations of spatial cooperativity. It has been found that, the 
model accounts for the presence of shear-rate fiuctuations in seemingly quiescent regions, 
and the local deviations from the macroscopic fiow curve when the inhomogeneity of the 
driving, quantified by the Babel number, is large. The comparison with experimental data 
also provided us with an estimate of the size of an elastoplastic block, N0 ^ 2, roughh 
in agreement with experimentally-measured values of the size of a shear transformation 
These successes of the mesoscopic model are encouraging for further studies of statistical 
aspects of fiow in complex fiuids. The description in terms of fiuidity diffusion is reasonable, 
but imperfect. 

The much larger deviations observed in the presence of rough walls, on the other hand, 
are not described by a model that imposes a fixed velocity at the walls. This points to 
the contribution of another physical mechanism to the deviations. We hypothesize that the 
bumps of the particles against the surface asperities due to wall slip, and the long-range 
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elastic deformations they induce, may be at its source. Further experimental and numerical 
studies of wall rheology will be needed to quantify this mechanism. 

AN & JLB thank K. Martens, D. Rodney, L. Bocquet and P. Chaudhuri for interesting 
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grant ERC-2011-ADG201 10209. 



[1] P. Shearer, Introduction to Seismology (Cambridge University Press, 1999) p. 272. 

[2] A. Argon, Acta metallurgica 27, 47 (1979). 

[3] P. Schall, D. A. Weitz, and F. Spaepen, Science (New York, N.Y.) 318, 1895 (2007). 

[4] H. Princen, Journal of Colloid and Interface Science 105, 150 (1985). 

[5] A. Amon, V. B. Nguyen, A. Bruand, J. Crassous, and E. Clement, Physical Review Letters 

108, 135502 (2012). 

[6] K. Nichol and M. van Hecke, Physical Review E 85, 061309 (2012). 

[7] M. L. Falk and J. S. Langer, Physical Review E 57, 7192 (1998). 

[8] A. Lemaitre and C. Caroh, Physical Review E 76, 036104 (2007). 

[9] J. Goyon, A. Colin, G. Ovarlez, A. Ajdari, and L. Bocquet, Nature 454, 84 (2008). 

[10] L. Bocquet, A. Colin, and A. Ajdari, Physical review letters 103, 036001 (2009). 

[11] P. Jop, V. Mansard, P. Chaudhuri, L. Bocquet, and A. Colin, Physical Review Letters 108, 
148301 (2012). 

[12] J. Goyon, A. Colin, and L. Bocquet, Soft Matter 6, 2668 (2010). 

[13] J. Baret, D. Vandembroucq, and S. Roux, Physical review letters 89, 195506 (2002). 

[14] G. Picard, A. Ajdari, F. Lequeux, and L. Bocquet, Physical Review E 71, 010501 (2005). 

[15] K. Martens, L. Bocquet, and J. Barrat, Soft Matter 8, 4197 (2012). 

[16] T. Mura, Micromechanics of Defects in Solids (Springer, 1987) p. 587. 

[17] R. R. Hartley and R. P. Behringer, Nature 421, 928 (2003). 

[18] Videos are provided as Supplemental Material. 

[19] We have discarded the curves corresponding to the lowest two applied pressures in |ll^, because 

we are not sure of the accuracy of the measurements of the lowest shear rate fluctuations. 

[20] P. Chaudhuri, V. Mansard, A. Colin, and L. Bocquet, Physical Review Letters 109, 036001 

10 



(2012). 

[21] L. Isa, R. Besseling, A. Morozov, and W. Poon, Physical review letters 102, 58302 (2009). 
[22] J. R. Seth, C. Locatelli-Champagne, F. Monti, R. T. Bonnecaze, and M. Cloitre, Soft Matter 
8, 140 (2012). 

Supplementary material 

FICTITIOUS PLASTIC EVENTS AS SOURCES OF MECHANICAL NOISE 

To mimick the impact of the "bumps" into wall asperities induced by wall slip, we intro- 
duce fictitious plastic events along the wall. More precisely, we now model both walls as 
lines of blocks and select a fraction of these blocks at random as mechanical noise sources, 
that is to say, they shall release a constant plastic strain e^^y^^^ by unit time. Figl5] shows 
the resulting local flow curves for a fraction of flxed mechanical noise sources of 1/3 picked 
at random, and ell^^^^ = ±4.5. 

VIDEOS SHOWING THE SPATIOTEMPORAL DYNAMICS 

Files Video Op 7 5. mpg and VideoOp^S.mpg are videos representing the succession of plastic 
events in a channel flow. 

The system flows along the horizontal direction, towards the right, and the flow is conflned 
between two horizontal walls, at the top and at the bottom. A plastic event is depicted as 
a coloured square that gradually fades into black as the plastic event comes to an end. 

The color of the plastic event represents its principal direction, ie the yielding angle: 
Blocks that yield in the macroscopic shear direction, ie e^^^ = and e^y ^ , are painted in 
red, while plastic events with e^^^ 7^ and t^^y = are represented in blue. A color gradient 
is applied between these limiting cases. 

Video Op 7 5. mpg shows a system of A^^ = 7 blocks in the crosswise direction and A^^^ = 64 
blocks in the streamwise direction (only half of the horizontal axis is shown). The stress 
imposed at the wall is = 0.75 in model units, and the other model parameters are identical 
to those used in the paper. The video covers a time span At = 33.5 in model units. 
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Figure 5. Local shear rate a (y) vs local shear rate 7 (y) (averaged on streamlines y = est) in the 
micro channel, when fictitious mechanical noise sources of intensity el^y^^^ = ±4.5 are added on a 
fraction (1/3) of blocks on the wah lines. cr^=(4) 0.36, (•) 0.48, (■)0.8, ( )1.0, (t)1.1 in model 
units. Solid line: macroscopic flow curve. 

Video0p48.mpg shows a system of A^^ = 16 blocks in the crosswise direction and A^^^ = 64 
blocks in the streamwise direction (only half of the horizontal axis is shown). The stress 
imposed at the wall is = 0.75 in model units, and the other model parameters are identical 
to those used in the paper. The video covers a time span At = 28.3 in model units. 
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